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> 

' Abstract. We present Infrared Space Observatory mid- and far-infrared spectroscopy of the merging galaxy 

1/^ , NGC 6240, an object presenting many aspects of importance for the role of star formation and AGN activity 

1/^ ■ in [ultra] luminous infrared galaxies. The mid-infrared spectrum shows starburst indicators in the form of low ex- 

' citation fine-structure line emission and aromatic 'PAH' features. A strong high excitation [O IV] line is observed 

, which most likely originates in the Narrow Line Region of an optically obscured AGN. NGC 6240 shows extremely 

' powerful emission in the pure rotational lines of molecular hydrogen. We argue that this emission is mainly due 
to shocks in its turbulent central gas component and its starburst superwind. The total shock cooling in infrared 

. emission lines accounts for ~0.6% of the bolometric luminosity, mainly through rotational H2 emission and the 

CL(' [O I] 63^m line. We analyse several ways of estimating the luminosities of the starburst and the AGN in NGC 6240 

JL I and suggest that the contributions to its bolometric luminosity are most likely in the range 50-75% starburst and 

H '. 25-50% AGN. 

, Key words. Galaxies: individual (NGC 6240) - Infrared: galaxies - Galaxies: active - Galaxies: starburst 



1. Introduction is among the sources with indications for an active galactic 

nucleus (AGN) in the mid-infrared emission line spectrum, 

The double nucleus merging galaxy NGC 6240 is a similar indications can be obtained from the low resolution 

relatively nearby (cz=7339km/s, d=98Mpc) mfrared- ^id-infrared diagnostic diagram of Laurent et al. (2000). 

lummous system. At = 6.05 x 10 Lq (see e.g. Sanders g^^^^^j ^-ray studies have firmly established the pres- 

& Mirabel 1996 for the defimtion of the 8-1000/im lu- ^^^^ of powerful AGN activity (e.g., Iwasawa & Comastri 

mmosity Ljr), it is almost an ultralummous infrared ^ggg. yignati et al. 1999), in fact occuring in both nuclei 

galaxy and thus allows detailed studies of many aspects (Komossa et al. 2003). Coexistence of star formation and 

of the phenomenon of (ultra)luminous infrared galaxies j^^^ ^^^-^-^^ ^ j^^^^^^ ^^^^ [ultra] luminous infrared 

([U]LIRGs). For the dusty system NGC 6240, the infrared galaxies. The proximity of NGC 6240 and the extensive set 

luminosity is identical to the bolometric luminosity to niulti-wavelength observations make it a unique target 

good approximation, the optical/near-mfrared and X-ray studying the quantitative contributions of the AGN 

regions observed to be minor contributors. In the m- ^^^^ formation activity in a [U]LIRG with a relatively 

frared selected sample studied spectroscopically with the powerful AGN. Such studies are not only of interest for 

Infrared Space Observatory ISO by Genzel et al. (1998) it understanding of the local [U]LIRG population, but 



also for the high redshift infrared and submm selected 
populations to which the [U]LIRGs likely form the closest 



Send offprint requests to: D. Lutz 

* Based on observations with ISO, an ESA project with in- , , , 

struments funded by ESA member states (especially the PI ° 

countries: France, Germany, the Netherlands, and the United A second important question relates to the evolution 
Kingdom) with the participation of ISAS and NASA. of merging galaxies like NGC 6240. Recent near-infrared 
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ObsJD 


Date 
UT 


AOT 


Proposal 


■^Ranges 


Duration 
h:mm:ss 


09700106 


22-FEB-1996 


S02 


RGENZELJvlPEXGAL2 


16 


5:06:22 


09700207 


22-FEB-1996 


S02 


RGENZELJVIPEXGAL2 


11 


2:20:14 


31101401 


23-SEP-1996 


S02 


ASTERNBEJEGALAXI 


2 


0:56:40 


31801036 


29-SEP-1996 


S02 


RGENZEL^PEXGAL3 


7 


2:53:52 


48101201 


ll-MAR-1997 


S02 


PVDWERF_SHOCKS 


7 


1:49:02 


83300305 


25-FEB-1998 


S02 


EEGAMI_H2 


3 


1:15:26 


83500308 


27-FEB-1998 


S02 


EEGAMI_H2 


3 


1:00:34 


84400107 


08- MAR- 1998 


S02 


EEGAMI_H2 


8 


2:06:38 


86600118 


30-MAR-1998 


S02 


DLEVINE_FILLFUP 


3 


0:42:02 


27801108 


21-AUG-1996 


LOl 


HSMITHJRBGALS 


n.a. 


1:09:04 


81000771 


02-FEB-1998 


L02 


HSMITHJRBGALS 


n.a. 


0:20:38 


81000872 


02-FEB-1998 


L02 


HSMITHJRBGALS 


n.a. 


0:20:38 


81000973 


02-FEB-1998 


L02 


HSMITHJRBGALS 


n.a. 


0:20:38 


81001074 


02-FEB-1998 


L02 


HSMITHJRBGALS 


n.a. 


0:20:38 


81600475 


08-FEB-1998 


L02 


HSMITHJRBGALS 


n.a. 


0:20:40 


81600676 


08-FEB-1998 


L02 


HSMITHJRBGALS 


n.a. 


0:20:40 



Table 1. Log of the ISO observations of NGC 6240 presented in this paper. The number of wavelength ranges includes 'serendip- 
itous' data obtained by a different detector in different order simultaneously with a requested range, but excludes data outside 
the nominal SWS wavelength bands. 



dynamical studies (Genzel et al. 2001, Tacconi et al. 2002) 
suggest that merging ultraluminous infrared galaxies as a 
class will evolve into moderate mass, disky-type ellipti- 
cal galaxies. NGC 6240 stands out in this group by its 
very high stellar velocity dispersion (Lester & Gaffney 
1994, Doyon et al. 1994) suggesting a very massive ob- 
ject. Spatially resolved dynamical studies using CO (cold 
gas; Tacconi et al. 1999), H2 (hot gas; van der Werf et 
al. 1993; Tecza et al. 2000), and stellar absorption fea- 
tures (Tecza et al. 2000) reveal complex, not fully relaxed 
dynamics with maximum velocity dispersion between the 
two nuclei. The extraordinary luminosity of NGC 6240 in 
the near-infrared rovibrational lines of molecular hydrogen 
(Joseph et al. 1984) may be partly linked to this complex 
dynamics, as well as to shocked emission in the superwind 
flow that is driven by the NGC 6240 starburst (van der 
Werf et al. 1993). 

This paper presents ISO mid- and far-infrared spec- 
troscopic observations of NGC 6240 and discusses impli- 
cations for its starburst and AGN activity, and the nature 
of the powerful H2 emission. 

2. Observations and Reduction 

We have retrieved from the ISO archive all observations 
of NGC 6240 that wore obtained with the ISO Short 
Wavelength Spectrometer SWS (de Graauw et al. 1996). 
Table 1 lists the basic data. All observations were obtained 
in the SWS02 mode which produces short spectral seg- 
ments at the nominal SWS resolving power (~2000), cen- 
tered on selected wavelengths. The data were reduced in 
the OSIA^ 3.0 software framework which adopts the SWS 

^ OSIA is a joint development of the SWS consortium. 
Contributing institutes are SRON, MPE, KUL and the ESA 
Astrophysics Division. 



calibration files of the offline processing software version 
10.1. Standard reduction procedures were used, putting 
particular emphasis on identifiying and eliminating indi- 
vidual noisy detectors and scans affected by signal jumps. 
Repeated observations of the same range were offset by 
small values to a consistent median flux. Single valued 
spectra were then produced by kappa-sigma clipping and 
averaging the many independent measurements inside a 
resolution element. A resolution element of A/2000 was 
used for this step. While leading to a slight degradation of 
the nominal SWS resolution, this is favourable for signal- 
to-noise reasons and does not affect the wide line pro- 
files of NGC 6240 noticeably. Figure 1 shows the segments 
containing the detected lines and some ranges providing 
important upper limits. The spectra were rebinned with 
R~2000 and are five times oversampled with respect to 
that resolution. Table 2 lists the derived fluxes (estimated 
to be accurate to 20-30%), limits, and line parameters. 

All observations were centered on the region of the 
NGC 6240 double nucleus (Fried & Schulz 1983). The 
SWS apertures (between 14" x 20" and 20" x 33") include 
both nuclei as well as the extended region of emission in 
the near-infrared rovibrational lines of molecular hydro- 
gen (e.g. van der Werf et al. 1993). Partial analyses of 
the SWS spectra of NGC 6240 using only a subset of the 
data with earlier calibrations and addressing only certain 
aspects are presented in Genzel et al. (1998) and Egami 
(1998; 1999). 

We have also reduced an ISO-LWS 43-196/xm full spec- 
trum (LOl mode) as well as 6 grating line scan observa- 
tions (L02 mode) of the [NIII] 57/um line (see Table 1). 
The data were first reduced in LIA^ version 8.1, based on 
offline processing software version 10.1 calibration files. 

^ LIA is a joint development of the ISO-LWS Instrument 
Team at Rutherford Appleton Laboratories (RAL, UK - the 
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Fig. 1. ISO-SWS spectra of the hnes detected in NGC 6240. Some of the ranges from which important upper limits are derived 
are also shown. 
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Fig. 2. Infrared spectral energy distribution of NGC 6240 from ISO observations. Continuous lines show the ISOCAM-CVF and 
LWS spectra, histograms the ISOPHOT-S spectrum. Small crosses indicate continua from the SWS spectra at > 16/im, where 
the relative zodiacal light contribution is small. Asterisks indicate the broad band photometry of Klaas et al. (1997,2001) 



In particular dark current subtraction and absolute re- 
sponsivity correction were done interactively. Following 
the LWS recipes for faint sources we have further reduced 
the data using ISAP-^ version 2.1. We deleted glitches and 
memory tails, and averaged multiple scans using ISAP's 
standard clip & mean option with bin sizes preserving the 
sample spacing. To produce the continuous full grating 
spectrum shown in Fig. 2 we have shifted the 10 LWSOl 
sub-spectra in offset mode by values between 5 and 20 
percent (30% for detector LW5), using detector LW04 (« 
140 to 170 ^m) as reference ("good") detector. With an 
average effective LWS beam size of 77" (Gry et al. 2002) 
all relevant regions of NGC 6240 are covered. Fluxes or 
upper limits are included in Table 2. All measured lines 
are unresolved at the resolution of LWS grating spectra 
(150 to 300, « 1500 km/s). 

We have reanalysed the ISOPHOT-S spectrum already 
presented in Genzel et al. (1998) and Rigopoulou et al. 
(1999) using PIA'' version 9.0.1. The spectrum was ob- 
tained on 1996 February 15 using rectangular chopped 
mode (chopper throw of 180"). The pure on-source inte- 



PI Institute) and the Infrared Processing and Analysis Center 
(IPAC/Cahech, USA). 

^ The ISO Spectral Analysis Package (ISAP) is a joint de- 
velopment by the LWS and SWS Instrument Teams and Data 
Centers. Contributing institutes are CESR, IAS, IPAC, MPE, 
RAL and SRON. 

* PI A is a joint development by the ESA Astrophysics divi- 
sion and the ISO-PHT consortium 



gration time was 512 s. Steps in the data reduction in- 
cluded: 1) deglitching on ramp level. 2) subdivision of 
ramps in four sections of 32 non destructive read-outs. 3) 
ramp fitting to derive signals. 4) masking of bad signals 
by eye-inspection. 5) kappa sigma and min/max clipping 
on remaining signal distribution. 6) determination of av- 
erage signal per chopper plateau. 7) masking or correction 
of bad plateaux by eye-inspection. 8) background subtrac- 
tion using all but the first four plateaux. 9) finally, flux 
calibration, using the signal dependent spectral response 
function. The absolute calibration is accurate to within 
20%. Figure 2 shows the infrared spectral energy distribu- 
tion of NGC 6240 combining the ISOPHOT-S spectrum, 
the ISOCAM-CVF data of Laurent et al. (2000 and priv. 
communication), the LWS data and published broadband 
photometry. 



3. Results 

The mid-infrared fine-structure emission line spectrum of 
NGC 6240 resembles the spectra of starbursting infrared 
galaxies (Thornley et al. 2000; Forster Schreiber et al. 
2001; Verma et al. 2003) in showing strong [Nell] 12.81/xm 
emission and a moderate excitation as probed by the ratio 
[Nelll] 15.55^m / [Nell] UMnm. It differs from atypical 
starburst spectrum in the enormous strength of the rota- 
tional emission lines of warm molecular hydrogen, and in 
the fairly strong emission in [OIV] 25.89/im, a line which 
is one of the strongest tracers of the mid-infrared spectra 
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Transition 


ARcst 


'^Obs 


Flux 


FWHM 




/im 


^rn 


10"^° 


km/s 


Ha 1-0 Q(3) 


2.4237 


2.4832 


1.55 


540 


Ha 1-0 Q(5) 


2.4548 


2.5159 


0.83 


604 


H Br/3 


2.6259 




<0.90 




Ha 1-0 0(3) 


2.8025 


2.8720 


1.59 


548 


Ha 2-1 0(3) 


2.9741 




<0.35 




Ha 3-2 0(3) 


3.1638 




<0.30 




Ha 1-0 0(5) 


3.2350 


3.3144 


0.56 


436 


H Bra 


4.0523 




<2.5 




Ha 0-0 S(ll) 


4.1813 




<2.6 




H+1,2,3+1 ^3,3 


4.350 




<1.9 




Ha 0-0 S(9) 


4.6946 


4.8114 


3.90 


600^ 


[Fc 11] 


5.3402 


5.4707 


5.46 


875 


Ha 0-0 S(7) 


5.5112 


5.6469 


7.27 


636 


Ha 0-0 S(5) 


6.9095 


7.0803 


12.20 


580 


H Pfa 


7.4599 




<0.6 




Ha 0-0 S(3) 


9.6649 


9.9021 


9.24 


641 


[S IV] 


10.5105 




<1.1 




Ha 0-0 S(2) 


12.2786 


12.5821 


4.35 


484 


[Nell] 


12.8136 


13.1288 


16.20 


630 


[NeV] 


14.3217 




<0.8 




[Nelll] 


15.5551 


15.9405 


6.31 


675 


H+5,0^4,3 


16.331 




<1.0 




Ha 0-0 S(l) 


17.0348 


17.4497 


7.83 


657 


[Fell] 


17.9360 




<1.6 




[S HI] 


18.7130 




<2.0 




[NoV] 


24.3175 




<0.6 




[S I] 


25.2490 




<1.2 




[OIV] 


25.8903 


26.5384 


2.33 


672 


[Fell] 


25.9883 


26.6322 


1.85 


653 


Ha 0-0 S(0) 


28.2188 




<2.2 




[S III] 


33.4810 


34.3091 


3.57 


503 


OH 1/2-3/2 


34.6164 




T < 0.5 


no absorption 


[Sill] 


34.8152 


35.6786 


26.0 


731 


[Fell] 


35.3487 




<2.4 




[Nc III] 


36.0135 




<2.6 




HaO 643 — »■ 53a 


40.3367 




<3.7 




[OIII] 


51.8145 




<14.0 




[NIII] 


57.3170 




<4.0 




[01] 


63.1837 


64.744 


72.00 


unresolved 


[OIII] 


88.3560 




<6.0 




[Nil] 


121.8976 


124.95 


2.3 


unresolved 


HaO 3i3 — > 2oa 


138.5272 




<1.7 




[01] 


145.5254 


148.99 


4.00 


unresolved 


CO 17-16 


153.2667 




<1.8 




[CII] 


157.7409 


161.598 


29.00 


unresolved 



Table 2. Line fluxes and limits measured with SWS and LWS for NGC6240. Upper limits for line fluxes are for an adopted 
FWHM of 600km/s (SWS) which corresponds to an unresolved line for LWS. ^ FWHM enforced for the gaussian fit used to 
derive the line flux. 



of AGN (Sturm et al. 2002), but very weak or absent in 
the spectra of starbursts (Lutz et al. 1998; Verma et al. 
2003). 

The overall spectral energy distribution and mid- 
infrared low resolution spectrum (Fig. 2) also resemble 
star forming galaxies, with prominent emission in the mid- 
infrared aromatic 'PAH' features. The rotational lines of 
molecular hydrogen are so strong that some of them are 
easily detected in the low resolution spectra. 



3.1. Detection of AGN Narrow Line Region emission 

The detection of [OIV] at ~15% of the [Nell] flux is 
clearly above the regime of starbursts and approaches 
values for AGN (Fig. 3). Under the assumption of sim- 
ple mixing of AGN and starburst emission, the AGN will 
then contribute little to [Nell] but will dominate [OIV]. 
A word of caution is in place hero about possible other 
origins of [O IV] . The very low level [O IV] emission seen 
in starbursts is not generally attributed to AGN but to ei- 
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Fig. 3. Fine-structure emission line ratios for NGC 6240, com- 
pared to a number of bright AGN from Sturm et al. (2002), 
the starbursts of Verma ct al. (2003, excluding low metallic- 
ity dwarfs), and the supernova remnants IC443 and RCW103 
(Oliva ot al. 1999a, 1999b). The ovorplottod grid represents 
simple mixing by a linear superposition of starburst, SNR, and 
AGN line spectra. The loci of the grid lines indicate where the 
contribution to [Nell] is 0%, 10%, 20%, ... for supernova rem- 
nants and for AGN. The percentage contribution to other lines 
will be different. Note that the four leftmost AGN themselves 
have evidence (PAH emission) for significant star formation 
inside the SWS beam. 



ther supernova or wind-related ionizing shocks or very hot 
stars (Lutz ct al. 1998, Schacrer & Stasinska 1999). Given 
the dynamics and the superwind activity of NGC 6240 
and the extremely strong and likely shock-related molec- 
ular hydrogen emission, it is conceivable that, compared 
to a normal starburst, ionizing shocks producing [OIV] 
are at a significantly elevated level as well. One way to 
probe for such a contribution is to look for other emis- 
sions that are strong in ionizing shocks. An obvious choice 
here is [Fell] which is amply emitted in partially ionized 
zones of shocks (e.g. Graham ct al. 1987; Hollenbach & 
McKee 1989) and additionally boosted by shock destruc- 
tion of grains (e.g. Jones et al. 1996; Oliva et al. 1999a, 
1999b) onto which Fe is normally depleted. Indeed both 
starbursts and supernova shocks observed with ISO show 
strong emission in the [Fell] 25.98/xm line that is neigh- 
bouring [O IV] . Similar ratios of the two lines are observed 
in starbursts and SNRs (Fig. 3). In AGN, however, the 
[Fell] emission from shocked or UV/X-ray irradiated par- 
tially ionized zones is (relatively) fainter compared to the 
strong [O IV] emission from the NLR. The [Fe II] emis- 
sion in NGC 6240 is fainter relative to [O IV] than in star- 
bursts or supernova remnants which suggests the AGN to 
dominate [OIV]. This is quantified in the simple mixing 
model shown in Fig. 3 which assumes the spectrum to 
be a linear superposition of a starburst, an AGN, and a 
SNR spectrum. At the location of NGC 6240, the contri- 
butions to [Nell] are then about 82% starburst, 14% SNR 
and 4% AGN. Because of the strongly different intrinsic 



[OIV] /[Nell] ratios, the contributions to [OIV] are very 
different: 7% starburst, 18% SNR, and a dominant 75% 
AGN. This nominal mix as well as the complexity of the 
excitation mechanisms for these lines clearly allow some 
of the strong [O IV] emission to be due to a high level of 
non-AGN shock activity, but we will proceed in the fol- 
lowing with the assumption that [O IV] is mainly tracing 
the AGN proper. 

The higher ionization potential NLR lines or coronal 
lines are weak or absent in NGC 6240. No near-infrared 
coronal line detections have been reported in the liter- 
ature. This coidd be due to extinction. A limit of <0.26 
can be derived for the extinction- insensitive ratio of [Ne V] 
24.3/im and the neighbouring [OIV] 25.89/im line. This 
places the NLR of NGC 6240 towards the low end but 
within the excitation range of the Seyferts studied by 
Sturm et al. (2002), who measured ratios in the range 
from 0.18 to 0.37. The data are consistent with [OIV] be- 
ing AGN dominated. 

Egami (1999) suggested an overall dominance of the 
AGN contribution to the NGC 6240 fine-structure line 
emission, not only for [OIV] but also for the low excita- 
tion lines like [Ne II] . To make this consistent with typical 
AGN NLR hue ratios (e.g., Sturm et al. 2002), he invoked 
a model where the ionization structure of the NLR fol- 
lows an onion-like layered pattern from high excitation in- 
side to low excitation outside, with obscuration selectively 
blocking the higher excitation inner parts. High obscura- 
tion specifically towards part of the NLR is conceivable, 
also given the large amounts of gas in the central region 
of NGC 6240. We consider an extreme effect of obscura- 
tion on the line ratios unlikely, however, because of the 
non-ideal structure of real NLRs (e.g. Capetti et al. 1997) 
which does not show perfect layering of ionisation stages 
as in idealized spherically symmetric photoionisation mod- 
els. In addition, such a scenario would either require the 
same peculiar geometry to apply for both AGN, or in- 
significance of one of them, which is in contrast to their 
modest hard X-ray flux ratio of ^3 (Komossa et al. 2003). 

The combined evidence favours a model where the 
lower ionisation lines are dominated by starburst emis- 
sion but [OIV] by the NLR(s) of the AGN(s). The ab- 
solute flux of [OIV] from the NLR may still be affected 
by obscuration. NLR obscuration cannot be constrained 
from the data that are lacking extinction-sensitive ratios 
which are purely tracing the NLR. In principle, the [O III] 
nondetections (Table 2) also might constrain properties 
of the obscured Narrow Line Region (e.g. [O III] 88/xm / 
[O IV] 26/im < 2.6). However, our ignorance about density 
and ionization parameter of this region prevents convert- 
ing this into an obscuration constraint. 

3.2. Properties of the Starburst 

At an excitation ratio [Nelll] 15.55/xm / [Nell] 12.81/im 
of 0.39, NGC 6240 is well within the range covered by nor- 
mal and high metallicity starburst galaxies (Thornley et 
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al. 2000; Verma et al. 2003). This is still true when lower- 
ing the ratio by correcting for the likely noticeable contam- 
ination of [Ne III] by AGN emission. Estimated from AGN 
[NeIII]/[0 IV] ratios (Sturm et al. 2002) and the observed 
[OIV] flux, about half of the [Nelll] 15.55/im flux could 
be due to a relatively low excitation NLR, whereas AGN 
contamination of the stronger [Ne II] line is much less sig- 
nificant. A rough estimate of the intrinsic starburst neon 
excitation ratio is hence closer to 0.2 rather than the ob- 
served total of 0.39. This is in line with the overall picture 
of an ageing starburst derived on the basis of near-infrared 
data, e.g. the low equivalent width of the near-infrared hy- 
drogen recombination lines (e.g., Tecza et al. 2000). 

In the gas-rich circumnuclear region of NGC 6240, ex- 
tinction is most likely significant. If possible, it should be 
estimated from mid-infrared tracers representative of the 
starburst activity probed by the ISO observations. The 
two most reliable estimators, recombination lines and the 
[SIII] lines are of limited use for the NGC 6240 data, 
however. Because of their low equivalent width, no re- 
combination lines are significantly detected in our data. 
Comparison of the [SIII] 18.7/um line (inside the silicate 
feature) with the 33.48/im line can provide another es- 
timate of obscuration which is effectively a lower limit 
based on the assumption of [SIII] being in the low den- 
sity limit. The 18.71/im line is indeed undetected, con- 
sistent with high obscuration, but the overall weakness 
of [SIII] in NGC 6240 (as in other high-metallicity star- 
bursts, see Verma et al. 2003) makes the limit on the ra- 
tio with [SIII] 33.48/im still consistent with the intrin- 
sic ratio, and thus prevents derivation of a lower limit 
on obscuration. Another way for obtaining an estimate 
for the obscuration is to compare the near-infrared Br7 
flux of Goldader et al. (1997), obtained in a 3" x 9" N- 
S aperture, with the [Nell] flux. The ratio of 0.012 is a 
factor of 2 lower than the one derived in matching aper- 
tures for Br7 and [Ne II] fluxes in M 82 (Forster-Schreiber 
et al. 2001). This suggests 10-15Av higher extinction to- 
wards NGC 6240 than the well-studied extinction towards 
M 82 (Ay ~5 for the screen case, Forster-Schreiber et al. 
2001), in total 15-20Av in the simplifying screen assump- 
tion. This corresponds to ~lmag extinction of [Nell] and 
[Nelll]. This estimate is very uncertain for several rea- 
sons. It is sensitive to uncertainties in the flux of the weak 
Br7 line that is difficult to measure accurately on top of 
the stellar continuum. Adopting the smaller Br7 fiux of 
1.15 X 10~^^Wcm~^ obtained by integrating over a 5" 
diameter aperture in the Br7 map of Tecza et al. (2000) 
increases Ay by another ft!9mag. The estimate also relies 
on the assumption of similar metallicities in NGC 6240 and 
M 82. Our obscuration estimate is consistent with an esti- 
mate from the beginning of this century of Nh w 2 x 10^^ 
(Ay sslO) for the atomic hydrogen column in front of the 
radio continuum sources (Beswick et al. 2001). Since this 
column rc;fers only to atomic gas, it is effectively a lower 
limit to which any molecular column would have to be 
added. 
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Fig. 4. H2 excitation diagram for NGC 6240, showing for the 
lines observed with SWS the upper level population divided 
by the level degeneracy, as a function of upper level energy. 
The level populations are derived from observed fluxes without 
extinction correction. 

The picture of a low excitation starburst at obscura- 
tion somewhat above or similar to that of M 82 is also 
consistent with the 4a detection of the low excitation 
[NII]122/im line, at a wavelength much less sensitive to 
obscuration. The ratio of 0.14 to [Nell] is slightly above 
the estimated ratio 0.12 in M82 (Colbert et al. 1999 
vs. aperture-corrected Forster-Schreiber et al. 2001). The 
small critical density of the [Nil] line could additionally 
affect this ratio. The limits on the higher excitation lines 
above 50//m ([NIII], [OIII]) are close to fluxes expected 
for an M 82-like excitation; slightly deeper spectra should 
detect these lines. 

3.3. Extremely strong emission of warm molecular 
hydrogen 

Even more than the near-infrared rovibrational transi- 
tions, the mid- infrared rotational transitions of H2 are 
extremely strong in NGC 6240. The observed rotational 
lines in the range 0-0 S(0) to S(ll) (Table 2) carry a to- 
tal flux of 4.5 X lO^^^Wcm^^. Adding an approximate 
1/3 for the unobserved among the pure rotational tran- 
sitions, this corresponds to a luminosity of the II2 rota- 
tional lines of 1.8 x 10^ L© which is ~0.3% of the infrared 
(8-1000/um) luminosity of the galaxy. A convenient way 
to visualize the properties of the H2 emission is an exci- 
tation diagram (Fig. 4, adopting transition probabilities 
of Turner et al. 1977, see also Wolniewicz et al. 1998). 
Thermally excited II2 at a single temperature and high 
density will show all transitions on a straight line in such 
a diagram. The concave curvature seen in Fig. 4 indicates 
a mix of temperatures, similar to excitation diagrams for 
starbursts and AGN observed with ISO (Rigopoulou et 
al. 2002). The rotational emission in NGC 6240 is, how- 
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ever, not only scaled up but also shows a stronger popula- 
tion of the higher rotational levels. At the low energy end, 
both the limit T>145K obtained from the limit on the 
S(1)/S(0) ratio and the excitation temperature of ~365K 
derived from the S(l) and S(2) lines are within the ranges 
derived from those lines in other galaxies observed with 
ISO. Comparing lines originating at higher and lower ener- 
gies, however, the flux ratio S(5)/S(l)=1.56 is larger than 
in any of the starbursts and AGN presented by Rigopoulou 
et al. (2002). 

The upper level of the S(3) line which falls near the 

center of the silicate feature appears somewhat under- 
populated in Fig. 4, suggesting significant extinction of 
Ag.e/ijm ~lniag towards the H2 emitting region. The lev- 
els in the first vibrationally excited state from which the 
1-0 lines arise appear underpopulated as well. Extinction 
may contribute here as well, since the lines from the first 
vibrationally excited state are at shorter wavelengths. For 
most commonly adopted extinction laws the obscuration 
at 2.4-3.2/im will be similar or less than the extinction 
at 9.6/xm, however. Less than or about one magnitude of 
extinction is not sufficient to explain the difference be- 
tween the population of the vibrational ground state and 
first vibrational excited state. Subthermal population of 
the vibrational levels at not too high densities may play a 
role. 

A total mass of 2 5 x 10^ Mq for the 'warm' (a few 
hundred K) molecular hydrogen can be estimated from the 
observed S(l) flux, and making the plausible assumption 
that the S(l) emission originates in gas that on average is 
50-100 K colder than the excitation temperature estimated 
from S(2) /S(l). This is a few percent of the total molecular 
gas mass estimated from CO observations (Solomon et al. 
1997, Tacconi et al. 1999), a fraction not uncommon in 
starbursts and AGN. More significantly, « 5 x 10^ Mq of 
gas at an excitation temperature of ~1100K are needed 
to produce the observed S(5) flux. 

Draine & Woods (1990) proposed a model where 
most of the near-infrared rovibrational H2 emission of 
NGC 6240 is due to molecular gas illuminated by shock- 
produced soft X-rays. Such gas would also contribute to a 
lesser degree to the observed rotational line emission. An 
observational test for this model is the expected emission 
of H^, e.g. at 4.35 and 16.33/im. We do not detect these 
lines but the derived 3(t limits (Table 2) are factors of 
19 and 4 respectively above the predictions. They do not 
constrain the X-ray illumination scenario, also considering 
the much lower predictions of Maloney et al. (1996). 

Radiative excitation of molecular hydrogen by UV il- 
lumination in moderate density PDRs and subsequent 
cascading leads to a strong population of the higher vi- 
brationally excited levels which will be found above the 
pure rotational lines of similar energy in an excitation di- 
agram (e.g.. Fig. 3 of Timmermann et al. 1996; Sternberg 
& Neufeld 1999, Fig. 6 of Draine & Bertoldi 1999). The 
hmits set on the 2-1 0(3) and 3-2 0(3) hues rule out a 
dominant contribution from such moderate density PDRs 
for the NGC 6240 H2 spectrum. The higher vibrational 



levels are clearly not overpopulated with respect to the 
lower ones. This agrees with conclusions that have been 
reached on the basis of near-infrared spectra (e.g. van der 
Werf et al. 1993, Sugai et al. 1997). 

Beyond excluding UV excitation in low density PDRs 
as the prime contributor, the thermal H2 spectrum ob- 
served in the SWS range holds little direct clue as to its 
sources of excitation. However, the detailed multiwave- 
length observations available for NGC 6240 make a strong 
contribution of shocks likely. The CO velocity fleld in 
the inner gas concentration is highly turbulent, suggest- 
ing rapid dissipation in this large gas mass (Tacconi et al. 
1999). Shocks that will occur under such conditions are ef- 
ficient sources of rotational H2 emission up to fairly high 
excitations, with the excitation diagram generally showing 
a smooth trend with level energy (e.g., OMC-1, Figs. 2 and 
5 of Rosenthal et al. 2000) . The [O I] 63/um line expected to 
be strong in shocks (HoUenbach & McKee 1989) is also ob- 
served to be prominent in NGC 6240 (Table 2). Its fiux is 
^2.5 times that of [CII] 157/xm, unusually high compared 
to other starbursts and IR-bright galaxies (e.g., Fischer 
et al. 1999, Colbert et al. 1999). This high [OI]/[CII] ra- 
tio is even above the relatively high [01] /[CII] ratios ob- 
served in some [C II] deficient galaxies (Malhotra et al. 
2001, Luhman et al. 2003) which are often ascribed to in- 
tensely heated (high Go/n) PDRs. The very high ratio for 
NGC 6240 is probably not more than partially of the same 
origin - at log{L[ciij/ L-pm)=-2.7 (sec also Stacey et al. 
1991, Luhman et al. 1998) NGC 6240 is not [C II]-deficient. 
Morphology, kinematics and excitation conditions of the 
near-infrared rovibrational II2 lines have been variously 
used to argue for shock excitation of these lines as well 
(e.g., van der Werf 1993; Sugai et al. 1997). Single veloc- 
ity/density C shock models (Kaufman & Neufeld 1996) 
may not be able to reproduce simultaneously all the rota- 
tional lines. This can be plausibly accounted for by a mix 
of speeds and by a contribution of PDRs to the lowest 
rotational transitions. 

Depending on shock conditions, cooling by the pure 
rotational lines of H2 will correspond to a few percent to 
tens of percent of the mechanical luminosity of the shock 
(e.g. Kaufman & Neufeld 1996). The mechanical luminos- 
ity corresponding to the observed rotational H2 cooling 
may thus be approaching IO^'^Lq. A refined estimate is 
possible considering other potential major shock cooling 
lines. We do not detect in the LWS spectrum lines of CO 
and II2O that can also be major coolants for shocks in 
molecular regions. Table 2 includes limits for CO and H2O 
lines in the most sensitive regions of the spectrum. These 
are factors ^5 less than the fluxes of the strongest rota- 
tional molecular hydrogen lines. We may use these limits 
and the ratio S(5) to 1-0 Q(3) to search for acceptable 
conditions in the range of shock speeds and densities sam- 
pled by the shock models of Kaufman & Neufeld (1996). 
These argue for, on average, relatively fast (V~30km/s) 
shocks in gas of preshock density ~10^ cm"*^. Under these 
conditions, the rotational lines of molecular hydrogen will 
radiate ^40% of the shock mechanical luminosity. 
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The strong and probably shock dominated [OI]63/Kni 
Une indicates that this picture is not complete. It radi- 
ates another 0.3% of the bolomctric luminosity but is not 
produced in significant quantities in the nondissociative 
(C-) shock models cited above, where oxygen is assumed 
mostly locked in molecules (but cf. Poglitsch ct al. 1996). 
Strong [O I] emission is however a key result of models of 
dissociative (J-) shocks (Hollenbach & McKee 1989). A 
single shock model is not adequate to reproduce a galaxy 
with a wide range of shock speeds and preshock densi- 
ties. Assuming that the H2 rotational lines and [01] are 
the strongest coolants of shocks in NGC 6240, the total 
mechanical luminosity of these shocks must be of the or- 
der 10^° L0 even if some fraction of those lines has other 
origins. 

If derived from the kinetic energy of the turbulent 
central gas concentration (~ 3 x 10^ M©, <j «150km/s, 
Tacconi et al. 1999), a shock mechanical luminosity of 
10^° L0 will dissipate this turbulent energy within a few 
million years, i.e. on the dynamical timescale of the cen- 
tral gas concentration. Additional contributions may come 
from the mechanical luminosity of the starburst 'super- 
wind' in NGC 6240, estimated to be of the order lO^^Lo 
(e.g., Heckman et al. 1990, Schulz et al. 1998). 

4. What are the contributions of star formation 
and AGN to the luminosity of NGC 6240? 

We discuss several approaches to quantify the luminos- 
ity of the starburst and the AGN(s) in NGC 6240. All 
these approaches directly relate a mid-infrared or X-ray 
observable to the infrared (or bolometric) luminosity. In 
this they differ from the purely mid-infrared diagnostic 
diagrams (Genzel et al. 1998; Laurent et al. 2000) where 
extrapolation from mid-infrared to bolometric has to be 
considered separately (e.g. Sect. 3.6 of Genzel ct al. 1998). 
The uncertainties will thus include the uncertainties of 
the extrapolation to bolometric luminosity. In addition 
to normal uncertainties in measTirements and methods, a 
consistent difficulty of several of these approaches is in- 
complete characterisation of obscuration. Effectively, this 
leads to a situation where a lower limit to the luminos- 
ity of a component can be set with some confidence, but 
where it is very difficult to assess whether more of this 
component is present at higher obscuration or whether 
a different component takes over. It is then more useful 
to quantify what is observed, effectively defining a lower 
limit, than to make possibly biased interpretations of the 
obscured components. It is also worth remembering that 
obscuration towards spatially distinct regions (e.g., star- 
burst, Narrow Line Region, AGN dust source, AGN X-ray 
source) is likely to differ considerably. 

Since modelling from first principles is currently un- 
reliable for several of the constraints discussed below, we 
follow an approach of quantitative but empirical compar- 
ison to starburst and AGN templates. For NGC 6240 as 
well as for the (dusty) comparison starbursts we assume 
the bolometric luminosity to be equal to the 8-1000/im in- 



frared luminosity, while the relation between these quanti- 
ties is explicitly included for AGN. For brevity, we do not 
include the radio emission of NGC 6240 in our discussion. 
The radio results summarized by Beswick et al. (2001) 
are consistent with coexistence of starburst and AGN, as 
outlined below. 

4.1. Star formation: Low excitation fine-structure lines 

With [Ne II] mainly excited by star formation, this line can 
be used to estimate the ionizing luminosity and bolomet- 
ric luminosity of the starburst. During the evolution of a 
starburst, the ratio of ionizing luminosity of the hot stars 
to bolometric luminosity will decay, as will the excitation 
observed in the mid- infrared fine-structure lines (see, e.g., 
the models presented in Thornley et al. 2000). When try- 
ing to estimate the bolometric luminosity of the starburst 
component from the observed low excitation fine struc- 
ture lines, it is hence important to choose templates with 
excitation similar to the ratio [Nelll] 15.55/im / [Nell] 
12.81^m «0.2 derived above for the NGC 6240 starburst. 
Of the starbursts observed by Verma et al. (2003) this 
ratio is between 0.1 and 0.4 for M82, NGC 3256, and 
NGC 3690/IC 694. For these galaxies, we compute a mean 
Zo5(X[Neii]/-t'iR) = —3.03. The values for the individual 
galaxies are between -2.89 and -3.22. To meaningfully com- 
pare with the IRAS FSC-based infrared luminosities, we 
have considered their spatial extent compared to the SWS 
aperture by adding the two apertures for the two pointings 
on NGC 3690/IC 694 which together fully cover the active 
regions of this interacting system (see Fig. 1 in Verma et 
al. 2003) and by doubling the observed flux for M 82 where 
the SWS aperture only partially covers the starburst re- 
gion (see Forstcr Schreiber et al. 2001). 

Using this empirical calibration, the observed [Nell] 
luminosity of 4.8 x 10^ Lq converts into a starburst bolo- 
mctric luminosity of 5 x 10'^"'^ Lq, i.e. the major part of the 
bolometric luminosity of NGC 6240. This would be an un- 
derestimate if the obscuration were higher than in the star- 
burst templates. While there is evidence for higher obscu- 
ration of the NGC 6240 starburst than in M 82 (Sect. 3.2), 
we stick with the conservative assumption of similar ob- 
scuration because of the uncertainty of this estimate and 
of obscuration estimates for the other templates NGC 3256 
and NGC 3690. The starburst luminosity would be over- 
estimated if the metallicity were much higher than for the 
templates, which are already at supersolar metallicities, 
however (Verma et al. 2003). It would also be overesti- 
mated if the AGN contribution to [Nelll] were less, mak- 
ing the starburst more highly excited with a correspond- 
ingly lower ratio of total to ionizing luminosity. 

4.2. Star formation: Aromatic emission features 

The mid-infrared aromatic 'PAH' emission features are 
observed over a wide range of interstellar medium con- 
ditions from the diffuse ISM to star forming regions but 
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Fig. 5. Low resolution ISOCAM 5-16^m spectrum of 
NGC 6240 compared to a scaled spectrum of M82 (top) and 
a schematic model consisting of the sum of a scaled M 82 
spectrum and an absorbed power law continuum (bottom). 
Obscuration of the continuum is needed for a satisfactory 
representation including the region of silicate feature. The 
schematic model is inadequate in the regions of strong emis- 
sion lines because of differences in the emission line spectra of 
NGC 6240 and M82. 

are not observed to be strong in the immediate vicinity 
of AGN. A number of studies have successfully decom- 
posed low resolution mid-infrared spectra of galaxies into 
three components: A 'PDR' component dominated by the 
aromatic emission features, an 'HIF very small grain con- 
tinuum steeply rising towards longer wavelengths and a 
hotter 'AGN' dust continuum (Laurent et al. 2000, see also 
Sturm et al. 2000, Tran et al. 2001, and Forster Schreiber 
et al. 2003). Dust obscuration can be important (Lutz et 
al. 1998a, Tran et al. 2001), with the additional compli- 
cation of ice features (Spoon et al. 2002). The NGC 6240 
spectrum (Fig. 5) is similar to starburst spectra which 
are a superposition of 'PDR' and 'HII', but there is an 



additional continuum extending down to 5/xm, indicative 
of an AGN continuum which we will discuss in Sect. 4.4. 
Subtracting such a continuum, the peak flux density in 
the 7.7/im aromatic feature is «0.45Jy. 

The physics and chemistry of the carriers of the aro- 
matic features and their transient heating by UV photons 
originating over a very wide range of stellar types is com- 
plex and not fully understood. This certainly provides for 
a less direct link from PAH luminosity to the bolometric 
luminosity of a starburst than from luminosity of the emis- 
sion lines originating in the medium ionized by OB stars 
to bolometric luminosity of the starburst. Nevertheless, 
it is worth trying to extrapolate from the PAH emission 
to the associated bolometric emission by using starbursts 
as templates. In order to minimize aperture effects in the 
comparison to the IRAS FSC-based infrared luminosity, 
we use objects for which PAH data are available cover- 
ing a major part of the IRAS beam and/or the starburst 
region^. Comparing the PAH 7.7/zm feature peak flux den- 
sity in Jy, and the 8-1000/im IR flux in Wm~^, we obtain 
a mean log(S7.7/FiR) of 11.84 with a dispersion of 0.25. 
Application of this relation to NGC 6240 results in a star- 
burst luminosity of 2 x lO^^L©. 

This is an underestimate if the PAH emission is more 
obscured than in the comparison starbursts. Indications 
for the corresponding relative weakening of the 8.6 and 
11.3/Ltm features (in the wings of the silicate absorption) 
are not strong, however (Fig. 5). More important for as- 
sessing the robustness of this estimate are variations of the 
aromatic feature emission with environment. It is known 
since the IRAS mission that the aromatic features de- 
crease when approaching a hot star or when going in- 
wards from PDRs to HII regions (e.g., Boulanger et al. 
1988; Verstraete et al. 1996). Correspondingly, trends are 
observed in 'normal star forming' galaxies between the 
average radiation field intensity, echoed in the 60/lOG/im 
flux ratio tracing the large grain temperature, and the ra- 
tio of PAH emission to total infrared luminosity (Dale et 
al. 2001). This introduces additional uncertainties when 
extrapolating from PAH features to total luminosity and 
may cause underestimates if much of the star formation is 
concentrated in a small region. 

4.3. AGN: High excitation lines from the Narrow Line 
Region 

Here, we use the Narrow Line Region [OIV] 25.89/im lu- 
minosity as an indicator of the AGN luminosity and use 
again an empirical calibration. To avoid biasing our result 
by star formation contributions to the bolometric lumi- 
nosity of comparison objects, we choose as comparison 
those among the Seyferts observed by Sturm et al. (2002) 
that have measured [O IV] lines and are not listed in their 

5 ISOCAM-CVF M82, M83, NGC 253, NGC 520, 
NGC1808, NGC3690, NGC4038/39, Arp236, NGC 7252 
(O. Laurent, priv. comm.) and NGC 3256 ISOPHOT-S 
(Rigopoulou et al. 1999) 
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Table 1 as having strong evidence for circumnuclear star 
formation. The mean log(L[oiv] /Lir) of these 16 objects 
is -3.01, with a dispersion of 0.29. Using this calibration, 
the observed [O IV] luminosity of 7 x 10 Lq converts to an 
AGN infrared luminosity of 7 x 10^*^ Lq. This is an extrap- 
olation based on the relation of [O IV] and infrared lumi- 
nosities in the comparison Seyferts. In Seyfert galaxies, the 
10-1000/im infrared range will represent very roughly 1/2 
of the bolomctric luminosity (Spinoglio et al. 1995), with 
significant emission emerging at shorter wavelengths in the 
near- infrared, visual, and UV to X-ray ranges. The shorter 
wavelength emission in the integrated SEDs of Seyferts as 
in Spinoglio et al. (1995) will not always be pure AGN 
emission because of host contributions, but we assume 
it is dominated by the AGN. For an obscured AGN in 
NGC 6240, it is then plausible to assume that much of the 
NIR/optical/UV/soft-X part of the AGN bolometric lumi- 
nosity would be absorbed and reradiated in the mid- and 
far-infrared. We hence double the estimate of the AGN 
contribution to the infrared luminosity to 1.4 x 10^^ L©. 

This is an overestimate if the non-AGN shock contri- 
bution to [OIV] were large in NGC 6240, contrary to the 
estimate in Sect. 3.1. While we have tried to exclude ob- 
jects with PAHs indicating strong star formation from the 
comparison sample, it may still be an overestimate if the 
IR emission of the comparison objects still has significant 
contributions from star formation. This is possible because 
of aperture effects between the PAH measurements and 
the FIR continuum measurements in currently available 
data. An example for this problem is NGC 1068 with very 
weak PAH in the ISOPHOT beam, and considerable star 
formation outside this beam but inside the IRAS beam 
(Telesco & Decher 1988, Le Floc'h et al. 2001). Similar 
effects are hard to firmly exclude for other less well stud- 
ied large galaxies because of the mismatch between the 
IRAS beam and the ISOPHOT-S aperture used to search 
for PAHs indicating star formation. High spatial resolu- 
tion both in mid-infrared spectrophotometry and in the 
far-infrared would be needed to minimize these problems. 
Our AGN luminosity is an underestimate if the NLR of 
NGC 6240 is significantly obscured even at 26/xm. Several 
magnitudes of visual extinction needed to hide the NLR 
in the optical do not yet have a significant effect, however. 

4.4. AGN: Mid-Infrared continuum 

Mid-infrared spectra of AGN show a strong continuum 
due to warm dust in the vicinity of the AGN, either in 
the putative torus or on larger scales e.g. embedded in 
the Narrow Line Region. The spectrum of NGC 6240 has 
excess continuum over a pure starburst spectrum which 
may either be AGN related or an increased contribution of 
the steep 'HIT continuum. The excess continuum must be 
obscured in order not to violate the observed spectrum in 
the region of the silicate feature. Because of this obscura- 
tion and the steep rise of a starburst-related 'H IP excess, 
an excess AGN continuum is best quantified shortwards 



of the PAH complex, at ~5.9/xm (Laurent et al. 2000). 
Subtracting from the observed spectrum of NGC 6240 a 
scaled M 82 spectrum reproducing the 6.2 and 7.7/im PAH 
features, we estimate the excess AGN continuum to be 
~0.07Jy at 5.9Atm (Fig. 5). 

Again we use AGN without evidence for circumnu- 
clear star formation (no or only weak PAH features ac- 
cording to ISOPHOT-S observations) as templates. From 
40 such AGN with z<0.1 and observations in the ISO 
archive we estimate a mean relation between 5.9/im con- 
tinuum flux density in Jy, and the IR flux in Wm~^, of 
log(S5.g/FiR)=11.70 with a dispersion of 0.29. Application 
of this relation to NGC 6240 results in an AGN infrared 
luminosity of 4.1 x lO^^L©. Again, we have to consider 
that for the comparison Seyferts a noticeable part of the 
bolometric will emerge at short wavelengths that are likely 
absorbed and redistributed to the infrared in NGC 6240. 
As for the estimate based on the [O IV] line, we apply 
an additional factor 2 correction to account for the dif- 
ference between infrared and bolometric luminosity in the 
comparison Seyferts, arriving at an AGN bolometric lu- 
minosity of « 8 x 1O^"L0. 

This is an underestimate if the AGN dust continuum 
in NGC 6240 is heavily obscured. Relatively little obscura- 
tion is sufficient to produce the silicate absorption needed 
for a spectrum fit (e.g. Ay ~20 in Fig. 5), but high values 
are not excluded. Another potential bias which can cause 
an underestimate of the AGN luminosity is due to the 
choice of PAH-free comparison objects. If unification re- 
lated orientation effects cause variations in the AGN mid- 
IR brightness (Clavel et al. 2000), then our comparison 
objects selected to have a bright continuum standing out 
of the host galax;y PAHs might be above average in mid- 
IR brightness compared to an orientation-unbiased sam- 
ple. In the absence of observations that can spatially sep- 
arate host and AGN at mid and far-infrared wavelengths 
this bias is hard to break. Again, we will have overesti- 
mated the AGN luminosity if the IR luminosities of the 
comparison objects still contain significant starburst con- 
tributions, outside the ISOPHOT-S beam but inside the 
IRAS beam. We will also have overestimated the AGN 
luminosity if some of the excess 5.9/im continuum is still 
HII region related. 

4.5. AGN: X-ray emission 

The X-ray emission of NGC 6240 holds important clues to 
the nature of its AGN. The soft X-ray emission (0.5-2 keV) 
is still dominated by extended thermal emission tracing 
the starburst superwind (Komossa et al. 1998; Iwasawa 
& Comastri 1998). The hard X-rays (2-10 keV) are domi- 
nated by AGN emission reflected by cold and warm matter 
(Iwasawa & Comastri 1998) from two AGN located in the 
two nuclei (Komossa et al. 2003). Above lOkeV, absorbed 
direct AGN emission has been observed with BeppoSAX 
and RXTE (Vignati et al. 1999; Ikebe et al. 2000). The 
question is to estimate the bolometric luminosity of the 
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Constraint 


Starburst 


AGN 


Low excitation mid-IR lines 


50-100% 




PAH features 


>33% 




High excitation mid-IR hnes 




>24% 


Mid-IR dust continuum 




>13% 


Hard X-ray emission 




10-100% 



Table 3. Constraints on the contribution of starburst and 
AGN activity to the infrared luminosity of NGC 6240. 



AGN from the observed direct and/or reflected AGN X- 
rays. 

Spectral fits to the complete X-ray spectrum of 
NGC 6240 suggest an intrinsic (corrected for absorption) 
2-10 kcV luminosity of 2.2 x lO^^ ^ 2.7 x lO^^ergs"! 
(Vignati ct al. 1999, Ikebe et al. 2000; corrected to the 
distance assumed in this paper). The main uncertainties 
involved in the analysis of the 2-lOkcV emission are the 
geometry of the X-ray reflector, and the potential effect 
of scattering by the absorbing matter. Similar uncertain- 
tics arise when extrapolating from the directly observed 
very hard (>10keV) X-ray emission, due to the uncertain 
spectral index towards softer X-rays. In a next step, the 2- 
10 kcV luminosity has to be extrapolated to the bolometric 
luminosity. Assuming the mean radio-quiet quasar SED 
of Elvis et al. (1994, their Fig. 10), L2_iokev/LBoi w0.09, 
with significant scatter around this mean relation. In the 
regime of Seyferts, L2_iokov/LiR is typically around 0.1 
for a major fraction of objects that are not absorbed in 
the 2-lOkeV range, with a tail to smaller ratios (Risaliti et 
al. 2001). We have verified this by comparing observed 2- 
lOkeV fluxes of z<0.1 Seyfert Is from George et al. (1998) 
and extinction-corrected 2-lOkeV fluxes of Seyfert 2s from 
Bassani et al. (1999) with infrared fluxes. We have re- 
stricted the samples to objects where the absence of PAHs 
in ISOPHOT spectra suggests small star formation contri- 
bution, and excluded Compton thick objects. The mean 
log(L2_iokev/LiR) for 18 objects (12 Syl, 6 Sy2) is -1.0 
with a dispersion of 0.29. Following our previous approach, 
we assume iiR « 0.5Z/bo1 for the Seyferts. 

In total, the AGN luminosity estimated from the X-ray 
spectroscopy is in the range 2 x 10^'' to 6 x lO^^ergs"^, 
or between 10% and several times the luminosity of 
NGC 6240, depending on the estimate of the X-ray lu- 
minosity and assumptions for the bolometric correction. 
The X-ray data suggest the AGN plays a significant and 
possibly dominant role but cannot pin down the AGN lu- 
minosity accurately. 

4.6. Combining the constraints 

The five constraints on the contributions of star formation 
and AGN in NGC 6240 discussed in the previous subsec- 
tions are summarized in Table 3. We assign highest weight 
to two constraints. The X-ray constraint suggests that the 
AGN is contributing a significant and possible dominant 
fraction of the luminosity, the latter situation being within 
the uncertainty of the estimate based on the observed X- 



ray emission rather than just being allowed as an unseen 
obscured component. The second key constraint is the 
one on the starburst luminosity derived from the low ex- 
citation fine-structure lines, which is based on relatively 
straightforward analysis of starburst H II regions. It is dif- 
ficult to envisage H II regions that are low excitation and 
still have a high ratio of ionizing luminosity to total lumi- 
nosity, which would be needed to make the starburst con- 
tribution to the bolometric luminosity insignificant. We 
suggest that the most likely range of contributions to the 
luminosity of NGC 6240 is 50 to 75% starburst and 25 
to 50% AGN. This estimate is based on direct compar- 
ison of mid-IR and X-ray observables to the bolometric 
luminosity. Reassuringly, the agreement with pure mid- 
IR diagnostics (Genzel et al. 1998; Laurent et al. 2000) is 
good. 

5. Conclusions 

Infrared spectroscopy of NGC 6240 detects signatures of 
powerful star formation activity and of an obscured AGN. 
We have used these results in conjunction with X-ray spec- 
troscopy to estimate the contributions to the total lumi- 
nosity of starburst (50 to 75%) and AGN activity (25 to 
50%) . NGC 6240 is a luminous source of rotational emis- 
sion of molecular hydrogen amounting to 0.3% of the bolo- 
metric luminosity, a similar amount is radiated in the [O I] 
63/Ltm line. Shocks due to the turbulent central velocity 
field and the superwind are likely to dominate these ex- 
traordinary levels of emission, and imply a mechanical lu- 
minosity of at least 10^° L©. 
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